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ABSTRACT 


The problem of determining the quantity of class oem 
laboratories and instructors to train sections of students 
attending numerous distinct courses in a school such as 
the Fleet Ballistic Missile School is considered. 

A procedure is developed for determining feasible 
schedules in order to graduate a fixed number of trainees 
over time while minimizing the cost of facilities mix 


required. 
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I. INTRODUCTION 


The problem of determining the capacity of a school το 
train sections of students attending ers dus came 
courses has been considered by Willingham [1] as an opti- 
mization problem. He used linear programming techniques 
for determining the maximum number of convenings of each 
foe Ol course taught atv the school which can take place 
during a year period subject to resource constraints and 
lower bounds on the number of convenings of each type of 
Course. 

ats problem arises wich the conversion Of EM ul 
marines from Polaris to Poseidon missile systems. There 
EM Che necessity of assuring a proper quantity and un. 
Of resources to instruct all required personnel at the Naval 
Guided Missile School. Here the primary resources of inter- 
est are laboratory facilities of various types with their 
associated equipment, classrooms and instructors of each 
specialty pertinent to Poseidon technical training [1]. 

The problem as stated from the point of view of those 
responsible for the funding and operation of the school is: 
"Given the planning estimates from BUPERS of the number of 
mersonnel who require taining Im each Туре ol ον ος ο υπ 
specified time period, determine what level and mix of re- 
sources is adequate to carry out the mission." The present 


work will consider this last aspect of the problem, 





attacking it under the assumption that care umi screen 
the convenings of each section will yield to a minimum 
requirement of resources for a given time period. 

The above idea can be expressed more precisely and in 
general as follows: An important function of any school 
manning is to estimate the number er тасу 55 пашезу 
laboratories, classrooms and instructors, that are needed 
Memsupports a specified program of IMsrructionzand Com ne у 
in general terms, the schedule to obtain the best employment 
of the resources. By a specified program we mean the type 
νι. τε: ο ρε taught, their curricular —Á а and 
the quantity of students assigned to each one. The problem 
E then, to find that schedule which wequires the least 
ost set of facilities. 

Since a partícular schedule dictates a required set 
ob facilities, it is necessary to have a measure of some 
emaracteristic in order to be able to compare it to other 
feasible schedules. A measure of effectiveness of a sched- 
Mee Will be defined as the cost involved in the installation 
DENT facilities associated with it. 

In general, specifying a measure of effectiveness in 
a scheduling problem is specifying a set of equivalence 
classes of schedules and a preference ordering among these 
Classes. For example, specifying the Cost Or Chee ae aS 
needed as the measure of effectiveness for a partir wiam 
problem means that: 


(a) A11 sehedules that have a cost C, are equivalent, so 


one is indifferent as to which schedule is selected. 


(b) A schedule with cost C. is preferred to a schedule 


3t 


witch cost С, Tenn duc I а C4 2° 


From the above discussion it can be seen then, thar 


ls strictly less thane 


the minimal cost scheduling problem is one of evaluating 
eemeaules and the associaved ο ο Ol ες ος ο μαμα 
mne search for the minimum m Tt is clear that tnereis 
ey large number of possible seheaules Tor tne problem. 
Bis also clear that most such sehedules are uninteresting 
for reasonable measures of performance and only a few are 
worth considering. 

Although the least eost facility mix scheaulinezprebulen 
easy to state and to visualize; no easy solution is 
available. Evaluating the set of all feasible solutions 
(represented by all possible eombinations of sehedules and 
facilities) is, for practical purposes, a gigantic quest. 
The objective of this paper is to illustrate a method of 
solving the problem applying the concepts developed in 
[2]. A procedure will be presented for approaching one 
ss lle of the equivalence class of schedules chat requires 
pue least cost of the resources for a given program of 
instruction for a given time period. 

A first attempt to obtain an acceptable solution to 
the scheduling problem can be made by evaluating the costs 
of random schedules which, consequently are random variables 
too, estimating the distribution of “the Cesvs. aie somo 


attempting to estimate the probability that, im 4 sanguemen 





size n, at least one will be in the коес л Р, OL 
the population. For a discrete гапаошуаг таре F m üI 
OF Wels casce em SS Jt Pr be the probability that a single 
observation has a cost value equal to r. Then the proba- 
bility that the smallest member of a sample of size n has 


OSC value less than risas given Be 
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where Р, 

For obvious reasons it is desirable to have an estimate 
of the lowest cost which is as close as possible to the 
true minimum eost of all possible sehedules. Thus Р, must 
mepresent a very small fraction of the lower tail of the 
distribution. Consequently the individual р;'5, near the 
extreme are very small thus the probability is small that 
the most extreme cases will be presented in any random 
sample of feasible size [4]. 

If a sufficient amount is known about the population 
distribution such that it is possible to calculate approxi- 
mately the probability that an observation is less than a 
specified value, not far from the extreme of the distribu- 
tion, then the number of further observations necessary to 
improve on this value could be estimated. However, before 
an improvement is obtained, a large number of Useless 


values may be found. "nese consiqerationsmim bra r sn 





rate at which useful information obtained by random sampling 
is much less when estimating extremes than when estimating 
means. 

As in any combinatorial problem the task of evaluating 
the distribution of the cost due to schedules is far irom 
being easy and is time consuming. The suggestion is, there- 
Rope, to employ some kind of Neuristic algorithm on the 
search for the least cost value. If a schedule is gener- 
ated at random and reductions in the cost involved are ob- 
tained by applications of the heuristic mechanism and the 
Process is repeated n times, then the n lowest costs gener- 
ated will tend to be a random sample of extreme values from 
when Samples originally generated from the parent distri- 
Клоп. Then a Weibull distribution can be fitted το the 
data. The Weibull has the characteristic of being indepen- 
sent of the parent distribution and has as one of Its own 
parameters the minimum value in question, [2], [3]. 

The approach to solve the problem presented in this 
study is first to develop a mathematical model to determine 
the relationships between schedules and facilities mix 
having the cost of the resources needed as the schedule's 
effectiveness measure. Thereafter, to generate a random 
starting date for each section of each course and to deter- 
mine the cost of this schedule as given by the model; then 
by means of an heuristic procedure applied on стар веште ние 
trying to reduce the cost involved ЕОС ЕОС и 


minimum cost schedule is found. Finally to use a Sample 





of such schedules to estimate the parameters of the Weibull 
distribution so that a decision can De made Tuba a 
stop or to generate additional random schedules and το 
repeat the estimation of the parameters. 

This thesis is organized in four sections. In Section 
II is developed a mathematical description of the behavior 
er a school which must train personnel. In such a3 scheer 
the courses taught can be of different lengths and each one 
сап be composed of one or more sections, depending upon the 
necessity of skilled personnel. It will not deal directly 
ПИЕ che problem of a time-table of meetings wmathin each 
week. It is assumed that each type of facility is available 
a fixed ce hours per week and the number of trainees 
required in a year for a given course (or specialty) is 
aivided into sections and they are scheduled to start at 
metime during the year, It is this scheduling problem 
vith which this thesis deals. 

Section III deals with elements of the theory of extreme 
memes and the Weibull distribution. This qiSscripucionsr 
the feature of being bounded and describes the behavior of 
the extreme values taken from samples independently drawn 
Brom any parent population. 

Іп Section IV is proposed a heuristic cost reducing 
procedure and a stopping decision rule, Even though the 
procedure is rather slow in finding a local minima and was 
not used in the examples discussed in Section V , it is 


presented here because, in је ашспов = ораши и cee / 


Ше 





the essential characteristics required by ап асы 
to be useful for this problem. 

In Section V the results of an example involving a 
hypothetical school are presented and the findings discussed. 
In Appendix An is presented i Lhe paqan S Os Lo 
that proposed in [1] that has the feature of being balanced, 
mno ncacensemomiacillituENSsSecNbDy the sections Of або еп 

саге throughout the complete"perioa. 

The computer program for the scheduling algorithm and 
brief description of it are included in Appendix Bs) ine 
program was run on the IBM 360/67 system in the "W. R. 
Church" Computer Center at the Naval Postgraduate School, 


Monterey, Galifornia. 


2] 
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The major assumptions in the model are: 


a. Course curricular requirements are given in the 


course syllabus. 


b. Each unit of a specified facility is available a 


fixed number of hours per week. 


сен The total number of Students requirineenstruct ten 
is nearly constant and small changes can be absorbed 


by varying the size of the sections. 


ae SAH section is a physical group of students receiving 
сито рес тлеа туре ot instruction. 


No restriction is placed in the order in which the 
courses are taught, i.e., nO course has preference 
to any other course. Any section of any course as 


equally likely to start at any time of the period. 


A. FORMULATION OF THE MODEL AND NOTATION 

MOD T = (i: 1:2....5.5h) Бе the set of all different 
Ies Of courses that are taught. Let N, be defined as 
ширег ог sections of course 1 which must. be taught 
during the period of one year. Let S, = (s: 1,2,...,М№,) 
ропе set of sections of course i, for all ie I. 

Define J = {j: 1l.2,... 010) 8:60] st orta NEDE 
paetacillities suech-as classrooms, 1аротасогтезша d оир 
SINS TFUELOrSs. Tn themodder ЕЕЕ re 


designate the classrooms. 


iz 





Let K = {k: 1,2,...,T} be the set of all working weeks 
in a specified period of time whiten as one year ante 
May be defined as the number of units of facility j required 
in week k. 

Let а; be the duration in weeks of course i and Er be 
ene total amount of facility required to teden course type 
w or all i c LI and j « J, j 7 l, Por example, Ass is the 
Amber Of lab—hnours of some laboratory or number of Teerures 
or eontact-hours required by the nn group ot instrucciones 
We assume that а. апа Е are predetermined and are given 
as exogenous parameters to the model by the syllabus of 
Salen course. 

sequencing requirements for individual topics or usage 
of laboratories within a determined type of course is not 
@ensidered, and it is assumed that the total requirements 
for facility j by course i is evenly used among the weeks 


that comprise the duration of the course. Then it is possi- 


ble to define: 


as the weekly demand of facility j by one section of course 
tor all ie T and ale en Ber 

Let D. be the weekly availability of one unit ot 
цасаа ау), Im hours, LOBA MA 5 Ls fixed sora 


КАРЕ К, 
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We would like to set up a Schedule such that all sec- 
tions of all courses are taught during the time period of 
interest and the costs involved are minimized. Since under 
the same conditions, the schedules should be the same from 
ame year to another, if a Section has nog finished in: 
unscpuctlon ac the end of the 5.59 week it will be included 
Бо пере ата та οἵ. τπε cycle for vie seam: nee 98 το οι 
пе instruction. 

Define Kis as the starting date of the n Section е 
course i, for all i e I and all s e N., and 1 < К, < Т. 
Note that there exist as many schedules as the number of 
r oLtilnet combinations of the Kis 95 since, for a given 
course the sections can always be relabeled the total number 


of different ways in which they can be scheduled is 


N4-1 
ҮЛ . Now for each section of course i there are WEN 
i ic 


T 








combinations among the sections of course i', for 3! = 
Ma .b, 1' # 1. Then the total number of vditrecens 


schedules is 


L 
T exp[ > Mas 
i-1 
Τ (1) 
ID XD 
i=l | 
Let ο be equal to 1 if sections Or οσοι πμ. 


Progress Gunnin? stig ig week and 0 otherwise. Then if 


κ - Kis < 0, the section has not starten сии я = 0; 


isk 


if O < Kk - Kis < daz, the section is in Novse a0 Weck soe 


14 





E Kis 2 di, the section has finished its 


iS Gi UC ος mud ee = 0. Note that having Ss sk defined in 


that way implies that e = ] 


κ зе о: ~ Sasktas 
when Kis - k. This insures that once a section has started 
its instruction it will continue week after week until it is 
finished. 


Summing the τ: а еке S; in any k produces the 


number of sections of course i being taught during week k: 


oe 


mor cut o ie cm (2) 
Then Tijg Mn represents the amount in hours per week that 
facility j is required by all sections of course i auring 
week k. Summing over all i e I will produce the total number 
of hours of facility j needed to accommodate all courses in 
week k, j Á 1. 

If a given facility is not to be used more than it is 


available in any week, then we have 


К Pit Bay < Mj DA, 3 x sa ЈЕ Апае 


hoe 


Where Mix is the number of units of j which must be provided 
am week κ. 

The number of classrooms required depends on the total 
number of hours of meetings other than in laboratories 
occurring in a week. Let G be the set of indices J not 


corresponding to laboratories and j #1. 


1 





L 
Then L 5 Pest Asy will represent the total number of 
јеб isl % 


classrooms-hours needed or the total number of meetings 
other than in laboratories occurring in week Kana mus pe 
less than ог equal’ to ИВА спе соса ааа аи 


lk 5 
classrooms-hours which must be provided in that week. 


L 
Σ το ο ο. D PME OMEN M MEN 
Then јеб 1=1 2] Ik 1k J 
Henceforth, M, = maximum (M gi K=l,2,...,T) represents 


the number of units of facility j required in the year in 
Order that the proposed schedule, defined by Kes к = ae s 


s = ly... Ny, HENO nb) por E 


LE COST STRUCTURE, AN OBJECTIVE FUNCTION 

Since in comparing two different schedules with the same 
memoer of sections and courses, the overall operational 
costs in a year will be essentially the same in both 
schedules, the objective will be to minimize the initial 
EOS. 

сонде ρε the set of starting deves бозиро ою 
all courses, which constitutes a schedule, thus € = 
E εν Π.Ε кате = ο ε 5.) . Given a particular 
senedule £ the facility cost can bezroundzrreme 


CE s= Νο РВЕ (3) 


iS 
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where C 


J 


If j is a classroom or laboratory e. 


as equipment, ои азер у oil ceo. mc cen 


is the initial cost of facility Jore и 


includes elements such 


In the present study 


the initial cost to obtain an instructor is а=ашпес нон = 


zero and therefore does not enter in the evaluation of a 


given schedule. 


A subjective cost could be attached to the 


Uunscructors if 1t was desired to study their influence in 


the model. ΜΙ 


Μ = 


and 


M. > Tiana 
J —À 


A defined 


ix 

The objective ís 
Ns minimum. It is 
paper. 
but only to estimate 


me) s ZR. 


maximum (M k: k = 


must satisfy the following ΤΕΤΡ 


Ju mcr ит 


< МОВА or all Jae J andike K 
EU е 

A Т 

integer, 


by equation (2). 


to find that schedule: monen ena 


this problem which we treat in this 


We do nos hope to find the minimum сос Сеи Кеа е у 
its value and find a schedule with 


This estimate will be obtained by using the 


properties of the extreme-values and the Weibull distribution. 


Í 





III. EXTREME VALUES AND THE WEIBULL DISTRIBUTION 


A. A USEFUL DISTRIBUTION 

In this section some properties of the extreme values 
and the connection to solving the minimum Coet го певана 
Problem will be discussed. Then, using the properties of 
mae Weibull distribution, a procedure Wil рерге зейге - 
estimate the minimum cost. 

A good presentation of the elements of the theory of 
extreme values is given by [5], here we only call attention 
to the fact that the lowest value in a random sample of 
size n, drawn from a parent population with cumulative = 


EM LPibution funetion F(') (given by A" - minimum (U, ,U 


in РІА 


T., U.) where the U, are sample values) is treated in [9] 
under the name Smallest Order Statistic. Weibull in [3] 


gives the cumulative distribution of T in ts most 


περα] form 


Р, (2) = 1 - ехр(-ф(2)) (4) 


min 


where the function $(z) must have the following 


[5,2]. 


characteristics 


i) Because in the minimal cost-scheduling problem there 
exists a true minimum cost such that the minimum valve 
obtained from the sample cannot бе белом τ, ΠΕ, 
minimum cost represents a lower bound for the extreme- 
value distribution. Let z be this true minimum value, 
chen 
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II 
O 


F (2) 
min 


AZ CC 


l -exp(-e(z)) ICE m 


ii) In addition #(z) should behave as B(Z in 


B is a parameter that depends upon the extreme-value 


- t), where 


spread. B must be greater than zero if BE τα, is a 


i 
distributions uneton: 


One distribution that satisfies the conditions stated 


above is the Weibull distribution (in the form given in 


ШО): 


Py (2) =1- exp(- [N (5) 
min 


where 

F(z) = The probability that Zin is less than or equal 
to a certain value 2, 

Z = The minimum cost criterion, 

min 

E = A constant equal to the minimum cost ot Che 
population, 

V = A Constant parameter indicating tche value of 
the variable such that the probability that Amin 
is equal to or less than v is approximately 
0.63 (referred to as the characteristic smallest 
value in extreme-value theory), 

k = A constant parameter indicating the shape of the 


distribution. The distribution оС о сара Еа 
skewed, Symmetrical or negacıveilvyssyercce 
depending on whether k is less than, equal to, 

OY Preaver taaan s o 


ky, 


There is no strong theoretical support иис Rp n 
distribution but empirical studies have shown "e merits 
of its applicability [2,3]. The function has the advantage 
of requiring no knowledge of the parent οποιο σου ους 
having as one of its own parameters, the boundary value of 
Te Pess 

As has been suggested before, the minimum cost obtained 
after applying a suitable heuristic mechanism to a random 
value drawn from the parent population can be regarded as 
an extreme-value obtained from this same population. If 
the same procedure is repeated n times a sample of n random ' 
and independent observations of the extreme-value, i.e. of 
Zin will have been obtained. It is then possible L O 
estimate the value of the lowest cost facility mix-schedule 
Pee iiteing the Weibull distribution to the data. 

The author is well aware that the minimum cost samples 
Ea is erete and trying to fit it to a continuous Чаа 
bution will be a source of error. However, for a moderately 
υ ρε sample, the approximation will be close enough so as 
to accept the estimated lowest cost as a lower bound for 
any extreme-value obtained. Аз а general rough guide only, 
ioe 30 3s a reasonable range of valves of a for applying 


approximations in most asymptovle roll. aa 


B. PARAMETER ESTIMATION OP THE WEIBULL DISTRIBUTION 


Let 2 А а 


in, min, >” min, be n independent and iden- 
tically distributed random variables. S Ze 


< 


1) 
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2 (2) SS 2 (y < Zen) be the observavions wer dered 
according to their increasing size from на 
minimum costs with cumulative distribution function as (5). 


Recall that the distribution of Zn should pe τς. 


in 
then what is possible is that several different combinations 
и startinewdatves, i.e. scheamees , we Nave ο ροής ес 


value. Then since observations can be present with differ- 


ing frequencies, define 


OD 
Ey A ol n 


where A Ck) πο ο Ποσο ιο ος hore. e (OO 


is the Corresponding 
frequency with which the nn value occurs, and m is the 
προς of different values that ETT takes. Then; accor Tinie 
to [7] an unbiased estimator of the cumulative distribution 
function of the Weibull distribution is given by: 

EG 


m 
Mec DEM. 
min (k) pez mtl 





F(k) = D 


and 


e) = n 


1 


им Я 


a 


Having the data in the proper form and the estimate of 
pne distribution function, there Forza ОС Усе 
task of estimating the parameters of the Weibull distribu- 
tion. Recall that this particular version one oe 


tion has three unknown parameters. Although there απο ο 


el 





accurate methods of estimating these parameters [2], we 
will present a method based on the double logarithm trans- 
formation of the distribution function (b Pa Bin e 
Eo. 1]: 


log(1og(1-F,) ^) с MM 


1 

This is the eguation of a straight line on Weibull 
probability paper. The slope of the line ís the shape 
parameter k. Although the parameter v is not self-evident 
from the graph it can be easily estimated as it will be 
shown later. This procedure or plotting can be repeated 
mer aliterent estimates of z, starting with ζ = 2⁄1) and 
lowering its value until the best straight line fit to the 
data is obtained by the least squares method. 

The derivation of the regression equations for esti- 
mating the Weibull parameters will be done for grouped 
data. Let 2 с) апа gU бедеге ж а Ж So e гор 
πμ 1,2, ... „mand the associated cumulative distribution 


be: 


“(ку m pli) 
5 - a (7) 





Por simplicity denime 


Y, = log (log (1-F‘*?))74) keo η (8) 
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l m3 
> 


The regression equation as given in [9] is: 


у = У + b(X.-X) (9) 


where I, is the estimate I, given A and b. The slope of 


the line is given by 


m m 
m i Ey ý T. 
k=1 k=1 
Σ X, Y, - 
= m 
Ν᾽ k=1 
m 
2 
a ( Σ X.) 
2 k=1 
> X, _ 
k=1 m 


Writing the regression line in slightly different form 


Y, = bX + (Y - bX), 


and comparing it with (4) we have b as the value of k, the 


estimator of the shape parameter. In addition, since 


а а а = none 


~ 


v, the estimator of the scale parameter is obtained from 


25 





ae exp(£ - X). 


Once the values of k and V have рееп орраштпе ока 
given t determine the sum of squares of the differences 
between the values obtained by (8) using (7) and the values 
estimated by the regression equation, (9), i.e. 


SQR - je 


Hu ota d 


Mes M 


k=1 


Repeat this step for different values of t such that 
g" < rt, until some value of t, say st yields the minimum 
ΙΙ squared differences about the yegression i GSE ne 
associated estimators of k and v, say k* and v* are the 
estimated parameters of the Weibull distribution and ¢* is 


SU m te of the lowest cost; ie τμ 5ππτεαπππτ 


EZ . . 
min 
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IV. A HEURISTIC ACCORTIEHN 


Having developed the model and examined the theory of 
сл геле Values, the remdininge steps an zero en 
heuristic cost reducing procedure and some kind of stopping 


decision rule. 


DO COST REDUCING PROCEDURE 

As stated before, once a random sehedule has been ob- 
used. s heuristic procedure designated Go find a schedule 
whose cost is a local minimum can be applied to it. The 
wollowing heuristie procedure is suggested for moving from 
a random schedule Er to another of lower cost, say Бү : 

l. For each type of facility there is 3 maxima 
Mer of units required by a given schedule. However, this 
maximum number of units is not required for every week of 
lames year. Arrange the facilities in order of increasing 
number of weeks during which the maximum number of units is 
required. The order in which they will be selected in the 
searching process is given by the above ordering. In the 
косе of а tie the higher priority should be assigned to 
What facility whose initial cost is higher. 

2. Beginning with the facility that has highest pri- 
erity, to determine which courses are related Со тапа 
order them by decreasing size of the amount of facility 


required. 
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З. The search continues by, successively к 
the course which requires the largest amount orm hi CE lc uM R 
say 1, and, with the exception of the . section, fix the 
Starting dates otuall seetionsr'ot all course ME DO MSN 
F r rcCing date for- section 1i, τες. let K, ç =k for all 
k = 1,2,...,T, and compute the cost of the moqifieq schedule; 
If for some value of k, this yields a cost which is less 
than the previous one, retain this as a possible starting 
He for section 5. 

4. Repeat the steps 1, 2, and 3 until no further 
improvement is obtained. 

ως Searchime procedure possible can be shortened 
се 16 125 very unlikely that a facility which is fully 
utilized more than fifty per cent or more can diminished 
Dy one unit. Consequently those facilities with usage 
ücctor greater than fifty per cent should not be taken inte 


ΠΠ ΠΠ the search procedure, neither зпошта those whoce 


maximum number of units is one. 


mee A STOPPING RULE 
By using the above procedure we obtain a set of sched- 


ules ES whose cost Zi 


η can be viewed as a random 


sample of extreme values. These are, then, used to esti- 
mate z = ee) the minimum cost value that could be 
achieved. When translated into a computer program it is 
necessary to incorporate a rule such that after ta nuni 
trials a decision can be made whether to Stop Or ШОЛ ος 


making more observations. 
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For instance, the following rule could beuc rE АЕ 
expected decrease in cost by improving the scheduling Is 
greater than or equal to the cost of one unit of the least 
expensive laboratory, make adđitional trials; otherwise 


"Op. Formally 


IS UE 
T min Zn 


2 Она 


Bontaimuezsearching. -Thernumberzef additional trials q, is 


given by 
η ΖΜΜ | 
саг min 
Pr 
ZMM 
f (ZMM-Z . (E))dF SIE os ron 
e min Lin j 


where ZMM is the lowest value obtained at the moment of 


decision, 
¿* is the estimated value of the minimum cost, 


dF 15 the probability density те С а не ΕΠ 


¿38 
min Weibull distribution, 


στ.) is the cost of one unit of the least expensive 
laboratory, and 


Co is the cost os one addc Toma NORRIS e pine ES 


computer cost, analyste ο ο 
It may still be in the designer's interest to examine 
additional decision parameters. Two of those are extensively 
explained by McRoberts [2]. Here we will give only some 


2 





illustration of the reasons why they are of interest. 
These two parameters are the cumulative density function 
E. and the value ZMM ~ z*. The value of the shape param- 
eter К, tends to describe the skewness οἱ the distributions 
as illustrated in figure l. Assuming that areas A and B are 
equal, the probability of improvement P, _ пау represent 
a very small absolute potential for a ЕИ k or пау repre- 
sent a large band of improvement for a bigger k respectively; 
then it may be desirable to identify regions ZMM - 27% and 
р | in the rejection range such that questionable zones 
M evaluated more closely. To do so MeRoberts [2] 
presents a plot of ZMM - cz* against the probability of 
improvement. 

Figure 2 from [2] presents such a plot which is based 
em considering an equicost curve as a function ol Fo 156 
2 S (ZMM-c*) SON By 
min 
Bensidering the equicost curve in isolation, the threshold 


CI(«) and C, and plotted in the F 


decision line may be represented as intersecting the C 
eurve at the value that may also be considered a threshold 


value of Zn - z*, the absolute range of improvement below 


in 
which further search will not be feasible [2]. 

Zone A. A cost value point falling in here leads to 
a decision for continuous searchin μας 


region is analogous to the statistical Type 
Ji errors 


¿one B. Above threshold value, continue searching. 
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k< 3.259 К 75.259 





C Ζπιιη(!) 2 mın(2) 


Figure 1. Effect of parameter k [2]. 
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Figure 2. Decision zones for continued search [2]. 





Zone C. This zone lies below the threshold level to 
the left of the feasible absolute range. 


Therefore stop search. 


Zone D. Here the probability of improvement may be 
low, but the absolute range of improvement 


is large. This implies continuing the search. 
Other decision parameters could be the Facility Utili- 
zation Factor (UFX) and the Usage Time factor. Both have 
been added to the computer program. Actually they are two 
measures of performance. The facility utilization factor 
UFX(-) is defined as the ratio between actual requirements 
in hours of facility (-) by all courses in the period in 
we stion and the total availability in hours of facility 
in the same period: 


N. 


Е 

ШЕ, -- πο. ιν 
νη 

ГОГА СЕЈ 


The usage time factor PMAX is defined as the number of 
weeks that the facility is used at its maximum rate divided 
by T = 50 weeks. 

Four cases are immediately self evident: 


a, UPA small, PMAX small; росе сос ПК оос 
possible to eliminate one 
unit of а ο πρ. 
consequently its cost by 
Searohwneg cud DC 
schedule. 
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oe UPX Tiree; 


САО Ета 


αυ lange, 


PMAX small: 


PMAX large: 


PMAX large: 


it is very unlikely nas 
case can: occur Te NL E 
not so large and PMAX very 
small, TL turns out Tomoe 


the previous case. 


that implies that it may be 
possible to make a better 

use of that facility ος 
menting the number of sections 
ol the Course Or Cour- Ccm Adr 
make use of I£. In the pre: 
sent conditions no improve- 
ments are obtained by further 
Search, 


implies good employment of 
that resource. No further 


improvement may be possible, 


E 





V, SOME NUMERICAL EXAMPLES AND DISCUSSION 


The search procedure presented in the previous section 
[unns out to be faimlyrtimepeosnsuming ornrtrihesecompurer а 
MESS ror research δρυ οσους r eE IL SINE ceSss any ше ате 
further investigations in that area as well as to refine 
bue stopping rule if it is possible. 

In order to demonstrate the method in Sections III and 
IV , some hypothetical problems were run by varying the 
mumber of Sections of each course, without using the seareh 
procedure. The random sample was obtained by picking the 
least cost schedule out of a sample of ten and repeating 
his step n = 100 times. 

The most economical schedule found was compared against 
am opposed schedule similar το that suggested Dy ла e nan] 
[1]. A description of that schedule is given in Appendix A. 

Table 2, below, shows the results of one of those exan- 
ples. It was first assumed that the availability of instruc- 
Bere 15 infinite, i.e. their initial cost was Zero. пат 
Корте was compared with the resultgssot the same Scheauiec 
assuming a subjective initial cost of assigning an instrue- 
Bor to the school, equal for all types oO 3mstpuccobDs MEN 
either case the proposed schedule was as good as or better 
than the most economical of the sample, even less than the 


estimated minimum cost inferred from the sample values. 
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Other runs were made for the same problem, this time 
varying the seeds for the random schedules ΕΙ emer 
obtaining the same results ο ρα οσο ος ee 

Table 3 shows the results obtained by varying the 
Меј OL Sections Ol, seach Course: 

[Geis interesting το point out char nozscheauttr ars 
found with a cost less than the one proposed in Appendix 
A for each case. This suggest the possibility of starting 
ша that schedule and then applying a search procedure 
seeking a minimum value. Of course this might not have 
been the case if the heuristic procedure had been used. In 
addition the sample size is extremely small compared to the 
parent population of all possible schedules as shown in (1). 

This thesis and the paper by Willingham both assumed 
that every section of every course used facilities evenly 
Ser tes duration. The method in this thesis does nov 
require that assumption to be made and could be easily 
modified to reflect more accurately the actual Тао И 
usages over time for each section. If this were done it 
is unlikely that a "balanced schedule" of the type described 
Mche Appendix could so easily be found. Thus the schedule 
proposed in the Appendix would appear less favorable 


compared to the method described here, 
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TABLE 1 


Common data for a Hypotheticalhloepeon 


course type 1 2 3 Ц 5 6 7 8 9 


duration D MEME C 9 5 ONDE 8 
Requirements Relative 
en facility Tr pa 
οσους ο 
Facilities 
| 9 
2 1 5915 ορ. 
3 5 6 5 2 175"0 
H Y E l 7 T 159070 
5 10 2 9 2 84.0 
6 5 19.0 
T e 55.0 
8 T (БО 
9 4 0 
10 23 18 1.50 
11 8 TO 
12 8 PO 
13 20 1:0 
14 18 l| О 
15 5 0 
16 e 10 
17 D 2 120 
18 T í 1.0 
19 H E T 
20 8 8 150 
21 3 3 I 
22 LORO 
20 JEN 


ποσο απ το ο ο ο ο το ο 
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ΠΕ к 
OU E S EPIS ai 2 3 H D 6 T 8 9 


NUMBER OF 
SECTIONS H о lo 6 20 8 8 5 


Problem 1 

Ineluding a subjective instructor cost (*): 
EESDOSESDESSCHEDUBEER (in Appendix A)T Cot O 

MOST ECONOMIC: (from sample) 1.222290 
STATISTICAL INFERENCE: Estimated minimum cost -1421.98. 


Weibull, shape parameter .60574. 


Problem 2 

KOC ine luding instructor cost: 

PROPOSED SCHEDULE: (by Appendix A) 1383.00 

MOST ECONOMIC: (from sample) 1383.00 
STATISTICAL INFERENCE: Estimated minimum cost 1379.86. 


Weibull shape parameter . 66359. 


(*) The differences between these two solutions was in 
che number of instructors required. rhe nunber ol 
units of each type of laboratory was the same. 
Since the problem is to minimize the resources 
required, the proper way to analyse the problem 
Ys to include a cost cf assienime anzinstmucter 
to the school. 
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TABLE 3 


Problem 3 


COURSE TYPE ih 2 3 ! J 6 Т 8 9 


NUMBER OF 

SECTIONS 8 Ц Ц 8 6 белн} 5) 9 
PROPOSED SCHEDULE: (by Appendix A) Cost 1039.00 
MOST ECONOMIC: (from Sample) 1042.00 


STATISTICAL INFERENCE: Estimated minimum cos UNE 


Weibull shape parameter .70306. 
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VI. CONCLUSION AND SUMMARY 


This study addressed the problem of obtaining a 
Schedule- in a training school for the minimum cost ο шо 
Ken litjes required, to suppert a given training plan. 
based on the quantity of students in each type of course 
taught. The problem was approached by making a mathematical 
Кече ов the seneduling-facilities interfacing, translated 
later into FORTRAN IV language. 

Making use of the theory of extreme-values, a 
sample of random minimum costs was obtained and this data 
was fitted to a Weibull distribution to estimate the 
Minimum cost that it is possible to achieve. 

Further investigation is required in order to 
come up with an efficient search procedure that assures a 
Kst and reliable algorithm such that a minimum cost sehequle 
is obtained by successive use of it on a random schedule 
drawn from the population of all feasible schedules. 

The model without the heuristic procedure was 
tested by running several programs varying the number of 
sections required in a period of a year. In each case a 
proposed schedule (see Appendix A) was better than the most 
economical schedule based on a sample obtained from random 
schedules using the extreme-value theory, although this 


result would probably not hold if the Courses eigee шы еш 
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assumed to use the facilities at a uniform rate over their 
duration. 

From the present results it appears that a good 
method would be to use a schedule of the type proposed in 
(1) and then perhaps to apply some improving procedure to 


that schedule. 
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APPENDIX A 


PROPOSING A SCHEDULE 


A similar schedule has been proposed in [1] with no 
i jm of being optimal bus wich The бло О рате ше 
facilities are used more or less evenly throughout the 
period under study. 

Assume the planning period and the number of sections 
of each course to be taught in that period (usually : year), 


MET БО weeks ала м5 for aii 1 e ijeere Known 


1 

Define the interval between convenings for the sections 
of the ps course by CS, - [T/N, ], largest integer. Since 
the right-hand-side is seldom an exact integer define 
DIF = T - N,CS, = 1. 

ак 

Assign to sections numbered 1 the starting date Kis = 1 

For all 1 e I. Then the (N, - 1) sections will be scheduled 


to start their training CS, weeks apart, except for DIF 


I 
sections that will be scheduled CS. + ] weeks apart; it 
does not matter which of the (N, = |) sections are selected: 
That will give a balanced schedule for each type of course 


in the sense of facility use. 
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APPENDIX B 


THE COMPUTER PROGRAM 


The program has been written in FORTRAN IV for the IBM 


360/67 system at the Naval Postgraduate School. However, 


the program is self-contained and can be used in any other 


machine using FORTRAN IV programming language. 


The program has been made in subroutines for easy 


ENC provements or changes. In general itrio Teon rocca 


of: 
(a) 


SUBROUTINE RMDATE which generates the set οἱ random 
starting dates aszindieated below. 

Because of the assumption that any section of any 
course can start its instruction at any week of the 
year, let the random variable Kis ре спе вата 
date of the gn secun ot courses nen tne 


Probability that Kes be equal to Keg = k is given 


by 


= = + 
O ауы 


Then the сера мека не та рони а о Ss 


"us 
РК. = Kig? E T 


e ‚| 


ει Ὁ if n s 
is – * 2 


otherwise. 
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As we can see, the distribution of Kis is the 
discrete version of the uniforme (ripa πο p 
to generate a random schedule it is only necessary to 


use the following modification to any uni ormie naom 


number generator [0,1]: 


Me [T - RN 4 0.5], largest integer, 


where 


Kis is the generated (random) starting date of 
t 


H ® * 
the s secvion ot course i, 


T is the length in working weeks of the period 
that has been considered; T = 50 weeks per 


year, 


RN is the uniform random number in the interval 
EC a 


SUBROUTINE COMPUT which computes, according to the 
model developed in Section I, the following values: 


. Number of units of seach ο» 55ο: 115; 
required to meet the specified program of 


Coning: 


2. Number of weeks in which the above are actually 


required. 


3. The total usage, in hours per year, of each 
facility "ружа Tresupses. 


4. The least cost schedule (ZMM) of the sample. 
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(c) The search procedure is incorporated to the main 
program by means ol several Ub ONU ME E 


l. SUBROUTINE SELECJ which arranges the facilities 
according to the procedure described и 
Section IV. 


2. SUBROUTINE SELECI which selects a course. 


3. SUBROUTINE REVAL which puts up-to-date the 
values obtained by SUBROUTINE COMPUT after a 


Change оп опе svartine бае ασ ος ο e 

(d) SUBROUTINE WEIBULL which computes the estimates 
indacacved an section Lil: namely: 

і. с location parameter or minimum cost estimate. 
2. k shape parameter. 
3. V Scale parameter. 

(e) SUBROUTINE RULE which decides whether to stop or to 
continue searching, This part of the program makes 
use of the IBM SUBROUTINE QG9 (integration of a given 
function by the Gaussian quadrature method, nine 
points formula) which has been added under the name 
ος ο Πο RULE: ~ This subroutine a ORe Terming 
how many additional trials should be done, if any. 

One of the features of the program is that one or more 
predetermined schedules can be tested in case there exist 
constraints in the scheduling period: 

The program also gives a write-up of each tested schedule 
and the most economical schedule, including characteristic 
values as utilization factors, usage percentage and weekly 


requirements of facilities in units ot СА ЕЕ ЕИ 
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All data is fed to the program by means NDA PAM M 
MENTS at the beginning of it, except οσο νο μοι e 
and the proposed schedules that are fed by means of Data 
Cards. The arrangement οἱ the дата δα зна e 


1. The first card has the value of two programas oni: ell 
parameters, ST and SNO. ST indicates the number of 
proposed schedules that the programmer or decision- 
maker wishes to test. SNO indicates whether search 


procedure is desired or not. 


2. The next set of cards are course requirements, one 
for each буре ог facility. 


3. Finally are the proposed schedule cards. Each card 


contains the section's starting dates of one course. 
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